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Abstract The effects of colour reconnection (CR) at 
e^e“ colliders are revisited, with focus on recently de¬ 
veloped CR models. The new models are compared 
with the LEP2 measurements for e+e“ —)■ W’'"W“ —)■ 
< 11 ^ 2 * 13^4 Eind found to lie within their limits. Prospects 
for constraints from new high-luminosity e“'"e“ colliders 
are discussed. The novel arena of CR in Higgs decays 
is introduced, and illustrated by shifts in angular corre¬ 
lations that would be used to set limits on a potential 
CP-odd admixture of the 125 GeV Higgs state. 

1 Introduction 

Multiparticle production in high-energy collisions often 
involves many contributing intermediate sub-sources. 
The cleanest such example is e+e“ —>■ W+W“ —> 
11 ^ 2 * 13^47 or its equivalent with a ( 7 */Z°)( 7 */Z°) inter¬ 
mediate state. A more tricky one is multiparton interac¬ 
tions (MPIs) in hadronic collisions, wherein a variable 
set of (semi)perturbative partonic collisions together 
with the beam remnants are at the origin of the subse¬ 
quent hadronization. 

In neither case can a first-principles QCD calcula¬ 
tion be carried out to describe the particle production 
process. Instead string or cluster models are used [1]. 
Both are based on an Nc —t oo limit [2], wherein each 
colour-anticolour pair is unique. Thus, in the string 
model, each quark is at the end of a string, whereas 
a gluon is attached to two string pieces and thus forms 
a kink on a longer string usually stretched between an 
endpoint quark and ditto antiquark [3]. In simple sys¬ 
tems like e“''e“ —> 7 */Z° —> qqg such principles give 
unique topologies, but for more complicated situations 
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ambiguities arise. When these can be associated with 
the presence of unexpected colour topologies we speak 
of colour reconnection (CR). The historical example in 
this spirit is the decay B+ = ub —>■ ncW*" —> (uc)(cs) —)■ 
(us)(cc) —)■ K+J/t/: —> [4], where we have 

used brackets in intermediate states to delineate sepa¬ 
rate colour singlet identities. 

Similarly, for e+e“ —W+W”, with W+ —> qiq 2 
and W“ —)■ q 3 q 4 , to first approximation the qiq 2 and 
I 3 I 4 systems hadronize separately from each other. 
Deviations from such a production picture could be 
parametrized as an admixture of alternative colour- 
reconnected qiq 4 and q 3 q 2 systems. Such CR was 
highly relevant in the context of the W mass measure¬ 
ment at LEP2 [5,6], where a potentially non-negligible 
uncertainty was predicted. This led to the development 
of dedicated studies aimed directly at measuring CR 
in hadronic W+W“ events [7-10]. The most extreme 
CR models could be ruled out, but not enough statis¬ 
tics was collected to definitely distinguish between the 
more moderate CR models and no CR [11]. Neverthe¬ 
less such moderate-model reconnection in about half of 
all events provided the best overall description. 

Modelling and testing of CR in hadronic collisions 
is rather more complicated [12,13]. And yet the case for 
it playing an important role is compelling, e.g. from the 
rise of the average transverse momentum with increas¬ 
ing charged multiplicity. Thus, given the predominance 
of hadronic colliders in recent years, first with the Teva- 
tron and now with the LHC, recent CR studies have 
rather aimed to address the more complicated issues 
arising there, and has led to the introduction of several 
new models [14,15]. These rely only on the distribu¬ 
tion of final state partons just prior to the hadroniza¬ 
tion, making them directly applicable also to e“''e“ col¬ 
liders. And even if the CR effects are expected to be 
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significantly smaller in e"'"e“ than in pp, this is com¬ 
pensated by a cleaner environment allowing for higher 
precision. On the one hand, it is therefore highly rele¬ 
vant to go back and check whether the newly developed 
models are consistent with the LEP2 data. On the other 
hand, it is useful to consider what further tests may 
come in the future. As an example, the recently sug¬ 
gested 100 km e+e“ collider [16] would produce 0(10®) 
W+W” pairs, resulting in a statistical uncertainty on 
the W mass below 1 MeV, e.g. from semileptonic de¬ 
cays e+e“ —> W+W“ —> qiq 2 ^r'^. With the calculated 
mass shifts in the original CR paper of the order 10-20 
MeV [5] as a reference, such a precision should make it 
possible to rule out many CR models, and also (hope¬ 
fully) definitely confirm the presence of CR effects. 

With the discovery of the Higgs boson [17,18], a 
new arena for CR studies opens up. The Higgs state 
is very narrow — the expected width is of the order 
of 4 MeV — meaning that it is very long-lived. There¬ 
fore hadronization of the rest of the event already hap¬ 
pened and the produced hadrons already spread out by 
the time the Higgs decays. That is, the Higgs itself de¬ 
cays essentially in a vacuum, and has no interactions 
with the rest of the event, be that in e“''e“ or pp col¬ 
lisions. Among its key decay channels we find W+W” 
and TPtP, however, and here history repeats itself: fully 
hadronic decays would be sensitive to CR between the 
two gauge-boson systems. The variables of interest here 
are not only masses but even more the angles between 
the four hadronic jets. Such angles can be modified by 
CR, a phenomenon which was noted e.g. in the context 
of top mass studies [14]. CR uncertainties thereby affect 
precision measurements of the Higgs properties, one of 
the primary purposes of future e+e“ colliders. To be 
specific, the SM predicts the Higgs to be a CF-even 
state, which is also observed to be strongly favoured 
compared with the CP-odd alternative [19,20]. Exten¬ 
sions of the SM Higgs sector, however, allows for the 
observed Higgs to be a mixture of both possibilities. 
One place to search for deviations from the predicted 
SM Higgs behaviour is precisely the angular correla¬ 
tions in hadronic W+W” (or TPtP) decays [21]. Hence 
CR could introduce a systematic uncertainty, and in 
this article we do a first study on the size of such un¬ 
certainties in various CR scenarios. 

This paper is organized as follows. The different CR 
models we will compare are briefly summarized in sec¬ 
tion 2. The three next sections contain studies on three 
different sets of observables, namely, the W mass mea¬ 
surement, section 3, the search for CR effects in W+W” 
events, section 4, and the Higgs CP measurements, sec¬ 
tion 5. The article ends with a few conclusions, section 
6 . 


2 The CR models 

Our current understanding of QCD does not provide a 
unique recipe for CR. Therefore the best we can do is 
contrast different plausible scenarios, and let data be 
the judge what works and what does not. In this ar¬ 
ticle we will compare four different CR models, which 
provide a reasonable spread of properties and predic¬ 
tions. Before briefly presenting each of these models it 
is useful to outline some of the basic issues that are 
involved. 

One key aspect is what role is given to colour 
algebra. To illustrate this, again consider e+e” —> 
w+w- —> qiq 2 q 3 q 4 - Erom the onset, qiq 2 form one 
singlet and q 3 q 4 another. In addition, there is a 1/9 
probability that qiq 4 and q 3 q 2 “accidentally” form sin¬ 
glets. In some models such accidental matches are a 
prerequisite to allow a CR. In this sense, these models 
are not really about reconnections but about a choice 
between already existing singlets. The alternative is to 
view CR as a dynamical process, wherein (infinitely) 
soft gluons can mediate any colour exchange required 
to form new singlets. The original non-accidental sin¬ 
glets define an initial state that actively needs to be 
perturbed to create alternative colour topologies. As so 
often, these two pictures may be viewed as extremes, 
and the “true” behaviour may well be in between, with 
a bit of each. 

Here another aspect enters, namely the role of geom¬ 
etry/causality. With a CT « 0.1 fm, the decays tend 
to be separated on a scale an order of magnitude below 
the typical hadronic size, the latter also being the size 
of the colour fields stretched between colour-connected 
partons. It would thereby seem that the W+ and W“ 
colour fields fully overlap, at least in the threshold re¬ 
gion where the W’s are not too strongly boosted apart. 
Introducing causality, however, the colour fields take 
some time to grow to full size (e.g. in the SK-I model 
described later). Meanwhile they drift apart, thereby 
only partly overlapping, and with an overlap that de¬ 
pends on the motion of all the string pieces from each 
W decay. In models where geometry is allowed to play 
a role there is also a natural decoupling of the two W 
decays at energies well above the threshold region, or if 
the W width could be sent to zero, and this should not 
be spoiled by the “accidental” singlets. 

Einally there is also a selection principle: if there 
are many potential reconnections in an event, which 
are the one(s) that actually occur? This could be at 
random or involve some bias. The most common bias is 
to make use of the A measure, which characterizes the 
total string length [22]. That is, the smaller the A, the 
better ordered are the partons along the strings. The 
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full A expression is rather messy, so a commonly used 
approximation is 


A = ^ In 

ij 



( 1 ) 


where the ij sum runs over all parton pairs connected 
by a string piece and mo is of the order of a typical 
hadronic mass. The average hadronic multiplicity of 
a string piece grows roughly logarithmically with its 
mass, so a reduction of A corresponds to a reduction of 
the “free energy” available for particle production. 

Among the four different CR models considered in 
this study, SK-I and SK-II were developed for W mass 
uncertainty studies at LEP2 [5]. The gluon move model, 
GM, was introduced as a simple model, among a few 
others, to study the effect of CR in top decays [14]. 
Finally, the QCD-based model, CS, was introduced to 
look for effects in soft QCD, especially baryon produc¬ 
tion [15]. The first two models are only applicable for 
the hadronic decays in diboson production, whereas the 
latter two could be used for any process. All of the mod¬ 
els are available in (recent versions of) Pythia 8 [23], 
the first two having been (re) implemented expressly for 
this study. That program also contains another CR 
model [13], used by default, that relies on the MPI 
structure of hadron collisions and therefore cannot be 
used in e+e“. All of the algorithms are applied after 
the hard primary process and the subsequent parton- 
shower evolution, but before the hadronization step. 
Typically this means that each W contains a handful 
of gluons, in addition to the primary qq pair, when CR 
is to be considered. 

Both the SK-I and SK-II model utilize the space- 
time picture of strings being stretched between the dif¬ 
ferent decay products of the two bosons. A reconnec¬ 
tion between two string pieces from different bosons 
is allowed only when these overlap in their space-time 
motion. Since such an overlap is assumed associated 
with the possibility for dynamical soft-gluon exchange 
between the two overlapping colour fields, there is no 
colour-factor suppression for reconnection. The two ap¬ 
proaches differ in their definition of what an overlap 
means, taking two extreme limits by analogy with Type 
I and Type II superconductors, which explains their 
names. In SK-I the strings are imagined as elongated 
bags, and the probability for a reconnection is pro¬ 
portional to the integrated space-time overlap between 
two string pieces. (Up to saturation effects to ensure 
that probabilities stay below unity.) This model con¬ 
tains one parameter that directly controls the overall 
strength of the CR, which made it convenient for ex¬ 
perimental LEP2 studies. For SK-II the string is consid¬ 
ered to contain a thin core, a vortex line, where all the 


topological information is stored, even if the full energy 
still is spread over a larger volume. A reconnection can 
only occur when the space-time motion makes two such 
cores cross each other. This model introduces no spe¬ 
cial parameters, and therefore gives unique predictions. 
(In both models one parameter is used to describe how 
the strings decay exponentially in proper time, and in 
SK-I additionally the string width is a parameter, but 
these parameters are almost completely fixed within the 
string model itself.) Normally only one reconnection is 
made, namely the one that happens first in proper time. 
By default this reconnection may either increase or de¬ 
crease the total A measure, but in the primed variants 
SK-T and SK-II' only reconnections that reduce A are 
considered. The SK models were tested at LEP2, where 
only the most extreme versions of SK-I were ruled out. 
For the SK-I model best agreement with data was ob¬ 
tained with parameter such that approximately 50% of 
all events contain a reconnection, as already mentioned. 


The gluon-move (GM) model was introduced to 
probe uncertainties in the top mass measurement, while 
still providing an overall good description of data. It is 
a very simple framework, in which the reduction of the 
A measure is at center, whereas neither colour algebra 
nor space-time geometry are considered at all. It con¬ 
tains two different types of CR, the gluon move one that 
gives the model its name, and a flip mechanism. In the 
former, the change in A measure is calculated if any of 
the gluons is moved from its current location between 
two colour-connected partners to instead be located on 
the string piece of any other colour-connected pair, fig. 
la. The move that lowers the total A measure the most 
is carried out, repeatedly until the minimum A measure 
is reached. The move step is quite restrictive, in that 
a string stretched between a q and a q endpoint will 
remain so; it is only the gluons in between that may 
change. Therefore an additional flip step is carried out 
after no more moves are possible. The flip mechanism 
flips the colour lines between two strings when this can 
reduce A, fig. lb, thereby mixing up also the string end¬ 
points with each other. (This is similar in character to 
what in another context is called colour swing [24].) A 
string is only allowed to do a single flip, to avoid the 
formation of gluon loops. The strength of the CR can 
be controlled by excluding a fraction of the gluons in 
the above scheme, or by requiring the A reduction in 
a potential move/flip to be above some minimal value. 
The parameters used in this study were tuned to de¬ 
scribe the LHC minimum bias data (although not quite 
as well as the default model). To allow more control, 
three alternative versions are considered in this article: 
only including the move mechanism, GM-I, only the 
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Fig. 1 Example of the gluon-move (a) and the gluon-flip (b) reconnections in the gluon-move model. The dashed lines represent 
the colour configuration of the partons. 


flip mechanism, GM-II, and the combination of both 
methods, GM-III. 


The S'C/(3)-based model, CS, is similar to the GM 
model, in that it also minimizes the A measure by 
doing flips between strings. But it differs in two ma¬ 
jor aspects. Firstly, it relies on the SU{3) colour rules 
from QCD, together with a space-time causality re¬ 
quirement, to determine whether two strings are al¬ 
lowed to reconnect or not. Secondly, it introduces a 
junction type of reconnection that is unique to this 
model. The use of SU{3) colour rules is a choice of 
philosophy, as already discussed. It limits which string 
pieces may flip with each other by requiring matching 
colour labels, i.e. that the colour flow is ambiguous al¬ 
ready by the colour assignments of the partons. It is 
possible to change the QGD-based default value, how¬ 
ever, in the extreme case such that all string pieces may 
flip with each other. For a flip between any two string 
pieces it is further required that they are in causal con¬ 
tact with each other, i.e. that each has had time to 
form before the other has had time to hadronize. The 
detailed formulation of this requirement is ambiguous, 
however, so a few options are available, with a tune¬ 
able parameter. The appearance of junction structures 
offers a clear extension relative to the other models. An 
(anti)junction is a point where strings stretched from 
three (anti)coloured quarks meet. In e+e” events they 
must be created in pairs, one junction and one anti¬ 
junction. When events hadronize, one (anti)baryon is 
created around each (anti)junction, thereby introduc¬ 
ing a new mechanism for baryon production. It is more 
important for high-energy hadronic collisions than it is 
for the studies in this article, however. A possibility not 
considered is that of colour ropes [25,26,24], where sev¬ 
eral parallel strings combine into one of a higher colour 
representation. If existing at all, ropes are more likely 
to play a non-negligible role in hadron or heavy-ion col¬ 
liders, where the beam axis offers a natural alignment 
of many strings. 


3 W mass measurements 

One of the key tasks of LEP2 was to determine the W 
mass, on its own right and as a test of the Standard 
Model consistency. Measurements were done both in 
the fully hadronic and in the semileptonic channels [10, 
27]. Both of them provide similar statistical errors, but 
the fully hadronic channel has a larger systematic un¬ 
certainty, due to the GR contribution. The uncertainty 
estimate depends on the analysis method as well as on 
the choice of CR models considered (and on their pa¬ 
rameters), but was found to be of the same magnitude 
as the statistical error. The large expected decrease 
in the statistical error at future e+e“ colliders would 
make the fully hadronic channel irrelevant for W mass 
measurements, unless the GR uncertainty could be con¬ 
strained by other means. This was already considered 
at LEP2 [10], where W mass measurements for differ¬ 
ent jet cuts were used to constrain the SK-I strength 
parameter. 

In this section we want to turn the table, and study 
how a precision measurement of the W mass difference 
between the fully hadronic and the semileptonic chan¬ 
nels would constrain CR models and parameter values. 
For this relative comparison a full optimization of both 
cuts and analysis methods is not required. Instead we 
will follow the method outlined in [5] to provide a sim¬ 
ple estimate of CR effects. 

To this end, one million e+e” —>■ W+W” —>■ 
ciiT 2 'l 3 T 4 events were simulated for each CR model. The 
events are required to have exactly four jets using the 
Durham jet algorithm [28], with a k± cut of 8 GeV. In 
addition the jets are also required to have an energy of 
at least 20 GeV each and be separated by an angle of 
0.5 radians. The four jets can be combined into two W 
bosons in three different ways. A few options for picking 
the “right” combination are considered: 

1. With the access to MG truth information, one can 

try to match each jet with a outgoing parton of the 

W decays. This is done by picking the match that 
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minimizes the product of the invariant masses be¬ 
tween each jet and its associated parton. 

2. One can use that the W mass is known to be close to 
80 GeV, and so minimize |mw ~ 80| to find the de¬ 
sired match, where m-yv is the average reconstructed 
W mass. 

3. Instead of requiring the average to be close to the 
known W mass, both masses individually could be 
optimized to be close to 80 GeV, i.e. minimize 

- 80| -b - 80|. 

4. At threshold the jets from the same W are almost 
back-to-back. A match can therefore be found by 
maximizing the sum of opening angles. 

To a large extent these methods pick the same com¬ 
binations, and thus they give similar results. Most of 
the problems arise in events with hard QCD radiation, 
where none of the methods are expected to work well. 

The W mass is calculated as the average of the two 
chosen W combinations. Since the target of this study 
is GR effects, the Breit-Wigner broadening of the mass 
spectrum is removed by subtracting the average of the 
produced W bosons. The results for all the methods are 
listed in table 1. The results for SK-I and SK-II differ 
slightly from the result in the original paper [5], which 
is due to the p_L-ordered shower in the newer versions of 
Pythia not being identical with the older mass-ordered 
ones of the time. 

The GM model shows an interesting behaviour; the 
move mechanism lowers the W mass, while the flip 
mechanism increases it, and the two effects accidentally 
cancel each other in the combined result. This may be 
understood as follows. If a gluon from Wi is radiated at 
a large angle, such that it will move closer to the decay 
products from W 2 , the move mechanism will connect 
the gluon to W 2 , fig. la. This will increase the mass of 
W 2 and decrease the mass Wi , but the decrease is larger 
than the increase, leading to the observed lower aver¬ 
age mass. The flip mechanism instead will connect jets 
between the two Ws, and thereby increase hadroniza- 
tion production of particles outside the W “cones”. This 
leads to larger opening angles, and thereby larger W 
masses. These two explanation will be revisited when 
studying the dedicated GR measurements. The com¬ 
plete cancellation is accidental, however, which becomes 
clear when the energy is varied. The SK-I and SK-II 
models also show opposite-sign effects, thereby further 
stressing the message that the mass-shift direction of 
GR effects cannot be taken for granted. Finally, the GS 
model shows no significant shifts, which will be a gen¬ 
eral trend throughout all the analyses. The limitation 
from the colour rules and the requirement of a lower A 
make effects very small at e“'"e“ colliders. By removing 
the colour constraints (GS max), the model starts to 


show an effect. This extreme case is already excluded 
at hadron colliders, however. 

A new collider should have the capacity to increase 
the energy beyond the W+W” threshold. And as was 
already observed for the SK-I model [5], the GR effects 
depend on the GM energy. There are two competing 
effects: firstly, the effect of a single reconnection be¬ 
comes larger with increased energy, and secondly, the 
probability to have two overlapping strings decreases 
with energy. The GR mass shifts for different CM en¬ 
ergies can be studied in tab. 2. Method 4 is here not 
included, since the maximum-angle method is only reli¬ 
able close to the threshold. The differences between the 
methods become smaller at higher energies, since the 
boost makes it easier to find the right combinations. 
The actual shifts increase at the intermediate energy, 
but drop when the energy is increased further. The only 
model that does not show this trend is the GS model, for 
which almost no effect is seen at any energy. The large 
shifts at the two higher energies for the other models 
provide a compelling argument to repeat the measure¬ 
ments at these energies. It should be recalled, however, 
that less statistics is expected at the higher energies. 


4 Four-jet angular distributions 

The direct searches for GR in W+W” events at LEP 
ruled out extreme parameter values for SK-I and po¬ 
tentially could also rule out some of the new GR mod¬ 
els. Especially the GM-I and GM-II models have that 
potential, since they were already observed to have a 
larger effect on the W mass measurement than the other 
models. 

The analysis relies on the particle multiplicities in 
the angular regions between two jets from the same W 
decay and from different W decays, respectively, to pro¬ 
vide a ratio that is sensitive to GR. The idea is that a 
reconnection will form a string between jets from dif¬ 
ferent W decays, thereby increasing the multiplicity be¬ 
tween those jets. In general, we will therefore expect 
the same-to-different ratio to become lower when GR 
is switched on. Several LEP experiments [8,9,7] per¬ 
formed this measurement. The results presented in the 
studies are after detector simulation, however, and as 
such are not directly comparable with the results ob¬ 
tained in this study. Instead we will rely on the ratio 
between the GR and the no-GR results (r, see later 
for exact definition), since detector effects are reduced 
for this observable. A preliminary combination of the 
different experiments gave r = 0.969 ± 0.011(stat.) ± 
0.009(syst.corr.) ± 0.006(syst.uncorr.) [29] correspond¬ 
ing to a 2.2 standard deviation disagreement with the 
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Method 

(Amy/) (MeV) 

{Smyf) (MeV) 

I 

n 

IT 

GM-I 

GM-II 

GM-III 

CS 

CS majc 

1 

-136 

-tl8 

-14 

-6 

-41 

+49 

+2 

+7 

+ 136 

2 

-73 

-tl3 

-13 

-7 

-28 

+34 

-1 

+3 

+73 

3 

-131 

-tl4 

-18 

-9 

-37 

+40 

-5 

+6 

+ 131 

4 

-1-131 

+ 10 

-18 

-9 

-27 

+31 

-3 

+3 

-131 


Table 1 Systematic mass shifts for the W mass at 170 GeV. The (/imw) value is the average reconstructed minus produced 
W mass for the no-CR baseline. The (Sfny/) is the additional shift for each CR model relative to this baseline. The statistical 
uncertainty on the latter quantity is 4 MeV. 


Method 

(<5mw> (MeV) {Ecm = 240 GeV) 

I 

n 

IT 

GM-I 

GM-II 

GM-III 

CS 

1 

2 

3 

+95 

+87 

+95 

+29 

+26 

+30 

+25 

+24 

+26 

-74 

-68 

-72 

+400 

+369 

+402 

+104 

+93 

+105 

+9 

+8 

+10 

Method 

{Smw) (MeV) {E^m = 350 GeV) 

I 

n 

IT 

GM-I 

GM-II 

GM-III 

CS 

1 

2 

3 

+72 

+70 

+71 

+18 

+18 

+18 

+16 

+15 

+16 

-50 

-50 

-50 

+369 

+369 

+369 

+60 

+60 

+60 

+4 

+4 

+3 


Table 2 Systematic W mass shifts at center-of-mass energies of 240 and 350 GeV, respectively. The (Smy^) is the mass shift 
in the CR models relative to the no-CR result. The statistical uncertainty is 5 MeV. 



Fig. 2 The (fir distribution at a center-of-mass energy of 183 
GeV. 


no-CR scenario. A later combined study [11] has in¬ 
creased this to disfavor the no-CR model at a 2.8 stan¬ 
dard deviation level, by combining with the mass shift 
results and performing a Ax^ fit. No separate r results 
were shown, however, and therefore we will have to rely 
on the preliminary combination. 

The event selection and analysis procedure varied 
slightly between the different LEP experiments. Two of 
the experiments relied purely on the angles to pair the 
jets [8,9], while one experiment also used the invariant 
masses [7]. We decided to mimic the analysis from the 
L3 collaboration [8]. A short recap of the event selection 
and analysis is presented, but for more details we refer 
to the experimental studies. 


The event selection requires each event to have ex¬ 
actly four jets with the Durham jet algorithm, with 
2/cut = 0.005. The two smallest of the six interjet an¬ 
gles are required to be below 100° and be non-adjacent. 
These are assumed to be the two regions between the 
different W decays, and are normally referred to as re¬ 
gions B and D. In addition two more angles are required 
to be between 100° and 140° and be non-adjacent. 
These are assumed to be the regions inside the W de¬ 
cays, and are normally referred to as region A and C. 
If several combinations are allowed, the one with the 
largest total opening angle is chosen. For each region 
the particles are projected onto the plane spanned by 
the two jets, and all particles are assigned a rescaled 
angle ipr = where is the angle from the parti¬ 

cle to one of the jets and (pjj is the angle between the 
two jets. This distribution is shown in fig. 2, where the 
different regions are separated by adding an integer to 
each. The final observable is defined as 

C"£-M + C)d^r 


Rn = 


-' 0.2 difr 
"0-8 dn 
0.2 dfpr 


{B + D)dipr 


The regions closest to the jets are excluded since they 
are mainly sensitive to the internal jet evolution. Fi¬ 
nally the ratio between the different CR models and 
the no-CR baseline is defined as r = Thus 

a deviation from unity would disfavour the no-CR sce¬ 
nario. The results for the various CR schemes are shown 
in tab. 3. As expected, all CR models, except for GM-I, 
predicts an r below unity. The GM-I model only allows 
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the gluon move reconnections, and therefore it does not 
reconnect the quarks at the string endpoints. Instead, it 
can take gluons emitted at large angles and move them 
to the other W string, thereby actually lowering the 
amount of radiation in region B and D, fig. la. This is 
the same explanation as for the lower W mass observed 
in section section 3. The GM-II model only does flips, 
which is exactly what this observable is optimized to 
measure. This is in fair agreement with observations, 
since this model shows relative large deviations from 
unity. The SK-I model with default strength is quite 
well in agreement with the actual measurement. For 
comparison the maximal SK-I model, where a recon¬ 
nection is always done, is also included. It gives too 
large shifts and so can be excluded. The SK-II mod¬ 
els and the CS model do not produce any large shifts. 
The maximal CS model, where the SU (3) rules are ig¬ 
nored and CR is only limited by the A measure, shows 
a larger effect and it can potentially be ruled out by ex¬ 
periments. It is, however, still relatively small compared 
to the other maximal CR models. 

In this study we consider several intervals, and not 
only the 0.2-0.8 considered in the original study. A 
clear trend shows that the smaller the interval, the more 
sensitive the observable becomes, i.e. varies more from 
unity. This is not surprising since the region closest to 
the jets are dominated by their perturbative behaviour. 
It should be noted that the statistics becomes worse for 
smaller intervals, but with the larger expected statistics 
at a new collider, a smaller interval than at LEP2 would 
most likely be preferable. 

To check if the new models are already excluded by 
the LEP measurements, the number of standard devi¬ 
ations from the measured result is calculated, tab. 4. 
The experimental uncertainties are assumed Caussian 
and added in quadrature. The only model excluded at 
the three a level is the CM-I model, which is the only 
model predicting a larger than unity r. The uncertainty 
is still too large to invalidate the other models, and a 
new collider with higher precision is needed to constrain 
these. 

The W mass measurement was seen to be more sen¬ 
sitive to CR at higher energies, and hence a similar 
effect is expected here. The method described above 
cannot directly be applied at higher energies, however, 
since the increased boost of the W bosons changes the 
angular distributions between the jets. Instead we ap¬ 
ply a method similar to method 3 in the W mass section 
to define the two angles within the W decays. The two 
other angles are defined to minimize the total sum of 
their angles. The results for the different energies are 
shown in table 5. The new method performs slightly 
worse at 183 CeV, i.e. the ratios lie closer to unity. 


This is especially evident when considering the maxi¬ 
mal CR models. At higher energies, however, the devi¬ 
ation from unity becomes larger for some of the more 
extreme models, indicating a better sensitivity, but this 
observable shows no sensitivity for the CS model. The 
moderate models do not show any significant variation 
with energy, and as such it is difficult to tell whether 
the potential limits on CR can be stronger at higher en¬ 
ergies. In general we expect a falling fraction of events 
with CR for higher energies, but more spectacular ef¬ 
fects for the events where CR occurs, so in the future 
we will need to search for more selective tests. 

As a slightly simpler observable, to test CR, it is 
possible to study the overall multiplicity. In most mod¬ 
els CR minimizes the A measure and therefore also low¬ 
ers the total multiplicity. This is normally compensated 
by a retuning of the hadronization parameters or the 
perturbative regime. But by comparing the multiplic¬ 
ity in fully hadronic and semileptonic W“*'W“ events, 
it is possible to directly probe CR. If no CR is switched 
on, the ratio 

is expected to be exactly equal to 2 (with £ = e or 
but excluding r). A simple study at a center-of-mass of 
170 CeV shows that indeed it is interesting to use this 
observable. Both the individual CM models show an 
effect, 1.96 and 1.97 for CM-I and CM-II, respectively. 
Contrary to the earlier observables, the two effects add 
coherently and the combined result is 1.93. With 1.97 
the CS model also shows more sensitivity in this observ¬ 
able as compared to the more complicated four-angle 
measurement. Similar results are also obtained for the 
SK models, so this would be an intriguing measurement 
for a future e+e” collider. 


5 Higgs parity measurements 

As discussed in the introduction, hadronic W+W” and 
Z°Z° decays of the 125 CeV Higgs offers a novel system 
for CR effects. Like in the W+W” studies above we 
should not expect big effects, so it is unlikely to be dis¬ 
cernible in the busy LHC environment. In a process like 
e+e- ^ 7 */Z 0 * ^ HOzO ^ E°e+£-, or /i+^- ^ H^ 
for that matter, detailed studies should become pos¬ 
sible, however, assuming sufficient luminosity. As be¬ 
fore, reconstructed masses and angles may become af¬ 
fected. Rather than simply repeating discussions along 
the lines of the previous two sections, we choose to il¬ 
lustrate possible effects for another set of observables, 
related to setting limits for CP violation in Higgs de¬ 
cays. We are aware that such tests can be performed in 
purely leptonic decays, say H —)• Z°Z° —)• ^+/r“e+e“, 
although with a much lower branching ratio. It can also 



inter¬ 

val 

pnoCR. 

_f_ 

I 

II 

IT 

GM-I 

GM-II 

GM-III 

CS 

I max 

CS max 

0.1-0.9 

1.1031 

0.9889 

0.9971 

0.9969 

1.0132 

0.9629 

0.9876 

0.9960 

0.9614 

0.9712 

0.2-0.8 

1.1482 

0.9802 

0.9916 

0.9931 

1.0293 

0.9440 

0.9918 

0.9910 

0.9360 

0.9781 

0.3-0.7 

1.1402 

0.9747 

0.9887 

0.9889 

1.0404 

0.9301 

0.9931 

0.9911 

0.9196 

0.9831 

0.4-0.6 

1.0883 

0.9702 

0.9823 

0.9880 

1.0460 

0.9181 

0.9882 

0.9920 

0.9068 

0.9810 


Table 3 Results for iJjv a-nd r for different intervals in ipr at a center-of-mass energy of 183 GeV. Two maximal CR models 
are included for SK-I and for the QCD based method, respectively. The statistical uncertainty on r is around 0.0025. 


1 1 no GR 1 I 

II 

IT 

GM-I 

GM-II 

GM-III 

GS 1 I max CS max 

1 1 2.0 1 0.7 

1.5 

1.6 

3.9 

1.6 

1.5 

1.4 1 2.1 1 0.6 


Table 4 Deviations from the measured result shown in number of standard deviations (rio- = (cexp — rth)/(^^’)exp). 


[GeV] 

pnoCR 

r 

I 

II 

IT 

GM-I 

GM-II 

GM-III 

CS 

I max 

CS max 

183 

1.9003 

0.9900 

0.9915 

0.9924 

1.0142 

1.0247 

0.9768 

0.9902 

0.9667 

1.0147 

240 

1.1764 

0.9820 

0.9935 

0.9933 

0.9857 

1.0130 

0.9362 

0.9993 

0.9030 

1.0006 

350 

1.4459 

0.9829 

0.9948 

0.9939 

0.9758 

1.0022 

0.9228 

1.0028 

0.8502 

0.9946 


Table 5 Results for i?jv and r for different center-of-mass energies for a fixed interval (0.2-0. 8). The statistical uncertainty 
on r is around 0.0015. 


be probed by the decay angles of the Z° produced in 
the association with the H° [30]. The purpose of this 
brief study is not to compare the relative merits of CP- 
violation tests in these different channels, but to stay 
with H —> W+W” —t qiq 2 q 3 q 4 and check what CR 
could mean there. To this end we will use a simplistic 
test on what could be the most sensitive variable. 
To simulate a mixed CP-even and CP-odd Higgs 
boson, we will use the Higgs doublet model already im¬ 
plemented in Pythia, with the option to allow CP- 
violation based on the expressions in [21]. We will as¬ 
sume that the 125 GeV Higgs is almost completely CP- 
even, with a small admixture of CP-odd. Allowing for 
an interference term between the two, the Higgs cross 
section can be written as 

fJ OC -f A^even^oddC , (3) 

where A,B,C depends on the kinematics of the event 
and the k determine the contributions to the different 
types. Since A, B and C are not of the same order of 
magnitude, a characterization in terms of a mixing an¬ 
gle is not convenient. Instead we use a definition based 
on the fraction, later referred to as parity fraction, of 
the events coming from either of the odd and the inter¬ 
ference parts of the cross section: 

r _ l^odd-^l + even fcoddCI 

^ " \kl,^A\ + \kl^^B\ + IftevenfcoddCI ’ 

For an almost CP-even Higgs, / = 0, this quantity 
provides a reasonable estimate of the amount of CP- 
violating interference introduced for the Higgs boson. 


The parity of the Higgs can be measured by study¬ 
ing the angles between the fermions from the boson 
decays. In the standard analyses of the spin/parity of 
the Higgs boson (see e.g. [31,32]), five such angles are 
defined, out of which three are sensitive to the parity 
of the Higgs. These three angles are: Oi, the polar angle 
of a fermion in the rest frame of its W mother, with 
respect to the direction of motion of the W in the H 
rest frame, 6 * 2 , similarly but for the other W, and <P, 
the angle between the two planes spanned by the decay 
products of the respective W bosons. The rest of this 
section will therefore be a study on the effect of CR on 
these three angles. 

To only have to consider the Higgs decay itself we 
have studied the process > H°, but this should 

only be viewed as a technical trick. All models are set 
up to easily handle this, whereas e^e“ —> H*^Z° would 
require a bit more bookkeeping for the SK models. Oth¬ 
erwise the models remain unchanged relative to previ¬ 
ous studies. The fact that at least one of the W’s have 
to be strongly off-shell implies that its lifetime is con¬ 
siderably reduced, and this is taken into account in the 
SK models. To estimate the effect of CR on the angles, 
100 million —)■ H° —> W+W“ —)■ qiq 2 q 3 q 4 events 

are simulated for each CR model and for each parity 
fraction, respectively. 

The events are required to have exactly four jets us¬ 
ing the Durham jet algorithm with a k± cut of 8 GeV, 
followed by an additional energy cut of at least 10 GeV 
per jet and a angular separation of 0.5. Two different 
methods to pair the jets were considered, either to max- 
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(c) 

Fig. 3 The three angles sensitive to the parity of the Higgs boson. Three different parity scenarios are shown together with a 
small selection of different CR models. 


imize the opening angles, or to minimize |Mw — 80| 
for a single W. The second method was found to be 
significantly more sensitive, and we will therefore re¬ 
strict ourselves to this method. The distribution for the 
three angles are shown in fig. 3. Deviations between 
the SM Higgs and the different parity fractions are vis¬ 
ible by eye for all the three angles. Both of the curves 
with nonvanishing CP-oddness show almost identical 
behaviours, indicating that the sign of the interference 
term is unimportant for these observables (at least for 
small deviations). Comparing the pattern of variation 
for the CP-violating models and the CR models, re¬ 
spectively, shows an interesting picture. For 9i and 02 
the deviations go in the same direction, whereas for <!> 
the deviations are in opposite directions. Thus a simul¬ 
taneous study in principle would allow to disentangle 
the two potential effects. 


To quantify the deviation from the no-CP-odd no- 
CR baseline, a simple test is applied to the distribu¬ 
tions. The most sensitive angle is di, and we therefore 
restrict our studies to this observable. A complete ex¬ 
perimental analysis most likely would combine all the 
angles in a multivariate analysis. For each parity frac¬ 
tion the is calculated, fig. 4. As expected the in¬ 
creases smoothly with this fraction. Similarly, the is 
also included for the different CR models. The crossover 
point is a simple indicator for when CR becomes an is¬ 
sue for Higgs parity measurements. This point occurs 
around 2-5%, with the higher values for somewhat more 
extreme CR models. Thus any limits significantly above 
this estimate can safely ignore CR effects. It should 
be stressed that also limits below 2% should be reach¬ 
able, once CR is carefully taken into account. This can 
involve (anti)correlations between the three angles, as 
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0 0.01 0.02 0.03 0.04 0.05 


parity fraction 

Fig. 4 Deviations between a CP-even Higgs without CR and 
models with either increased CP-oddness or a CR model. The 
deviation is quantified as the x^/NDF deviation for the (pi 
angle. 

already noted, but also studies of particle production 
patterns between the jets, like the one in section 4. 

6 Conclusions 

In this article we have studied the effects of CR at 
e+e“ colliders, with emphasis on fully hadronic W+W” 
events. We find that some newer models, implemented 
to study CR effects at hadron colliders, show different 
behaviours for e+e“. The CS model gives rise to very 
limited variations, whereas for the GM models one spe¬ 
cific scenario even shows large enough deviations to be 
excluded by the LEP data. 

Even if the concept of CR is quite straightforward, 
it allows for several different mechanisms to be at play. 
These potentially act in opposite directions, making 
interpretations difficult. This is clearly illustrated by 
the GM models, where GM-I predicts a smaller recon¬ 
structed W mass and GM-II a larger one. This high¬ 
lights the need for studying multiple models using sev¬ 
eral observables, to disentangle what is going on. Much 
further work is needed, but the outcome of the current 
simple study is fairly optimistic: given enough luminos¬ 
ity, at a few different energies, e“'"e“ should offer insights 
into CR mechanisms that complement those obtainable 
at hadron colliders. This complementarity between the 
“clean” e“'"e“ environment and the “dirty” pp one may 
hold the key to a deeper understanding of CR. 

The e+e” W+W” channel is not the only e+e” 
process where CR effects may be relevant. As an ex¬ 
ample we studied a Higgs parity measurement in the 
H —>■ W+W” —>■ qiq 2 T 3 T 4 channel. The variations from 
CR were of the same size as the introduction of 2-5% 
CP-oddness into the CP-even Higgs, depending on the 
choice of CR model. The main lesson is not the precise 


number for this particular observable, but to highlight 
the need to be aware of potential CR uncertainties for 
any nontrivial hadronic final state. 

Plans for future e+e“ collider usually include the 
possibility to reach the tt threshold. Then hadronic fi¬ 
nal states will start out with three colour singlets: one 
W from each top decay, plus one encompassing the b 
and b from the two decays. Like for the W+W“( 7 */Z°) 
background this increases the possibilities for CR ef¬ 
fects. Some early studies are found in [33], but updated 
and extended studies should be performed, including 
the new models. At the very least, it will be needed 
in order to estimate the expected CR uncertainty in 
the measurements of the top properties for possible fu¬ 
ture colliders. Many of the necessary tools are already 
in place in Pythia 8, although e.g. the administrative 
machinery in the SK models needs to be extended ap¬ 
propriately. 

Acknowledgments 

Work supported in part by the Swedish Research Coun¬ 
cil, contract number 621-2013-4287, and in part by the 
MCnetITN FP7 Marie Curie Initial Training Network, 
contract PITN-GA-2012-315877. Peter Skands is ac¬ 
knowledged for useful comments. 









11 


References 

1. A. Buckley, J. Butterworth, S. Gieseke, D. Grellscheid, 

S. Hoche, et al., “General-purpose event generators for 
LHC physics,” Phys.Rept. 504 (2011) 145-233, 

arXiv:1101.2599 [hep-ph]. 

2. G. ’t Hooft, “A Planar Diagram Theory for Strong 
Interactions,” Nucl.Phys. B72 (1974) 461. 

3. B. Andersson, G. Gustafson, G. Ingelman, and 

T. Sjostrand, “Parton Fragmentation and String 
Dynamics,” Phys.Rept. 97 (1983) 31-145. 

4. H. Fritzsch, “How to Discover the B Mesons,” Phys.Lett. 
B86 (1979) 343. 

5. T. Sjostrand and V. A. Khoze, “On Color 
rearrangement in hadronic W-|- W- events,” Z.Phys. 

C62 (1994) 281—310, arXiv:hep-ph/9310242 [hep-ph]. 

6. T. Sjostrand and V. A. Khoze, “Does the W mass 
reconstruction survive QCD effects?,” Phys.Rev.Lett. 72 
(1994) 28—31, arXiv:hep-ph/9310276 [hep-ph]. 

7. OPAL Collaboration, G. Abbiendi et al, “Colour 

reconnection in e-|- e-W-b W- at s**(l/2) = 

189-GeV - 209-GeV,” Eur.Phys.J. C45 (2006) 291-305, 
arXiv:hep-ex/0508062 [hep-ex]. 

8. L3 Collaboration, P. Achard et al., “Search for color 

reconnection effects in e+e“ —>• W~^W~ hadrons 

through particle flow studies at LEP,” Phys.Lett. B561 
(2003) 202—212, arXiv:hep-ex/0303042 [hep-ex]. 

9. DELPHI Collaboration, J. Abdallah et al., 
“Investigation of colour reconnection in WW events 
with the DELPHI detector at LEP-2,” Eur.Phys.J. C51 
(2007) 249—269, arXiv:0704.0597 [hep-ex]. 

10. ALEPH Collaboration, S. Schael et al., “Measurement 
of the W boson mass and width in e+e“ collisions at 
LEP,” Eur.Phys.J. C47 (2006) 309-335, 

arXiv:hep-ex/0605011 [hep-ex]. 

11. ALEPH, DELPHI, L3, OPAL, LEP Electroweak 
Collaboration, S. Schael et al, “Electroweak 
Measurements in Electron-Positron Collisions at 
W-Boson-Pair Energies at LEP,” Phys.Rept. 532 (2013) 
119—244, arXiv:1302.3415 [hep-ex]. 

12. T. Sjostrand and M. van Zijl, “A Multiple Interaction 
Model for the Event Structure in Hadron Collisions,” 
Phys.Rev. D36 (1987) 2019. 

13. T. Sjostrand and P. Z. Skands, “Multiple interactions 
and the structure of beam remnants,” JHEP 0403 
(2004) 053, arXiv:hep-ph/0402078 [hep-ph]. 

14. S. Argyropoulos and T. Sjostrand, “Effects of color 
reconnection on tt final states at the LHC,” JHEP 1411 
(2014) 043, arXiv:1407.6653 [hep-ph]. 

15. J. R. Christiansen and P. Z. Skands, “String Formation 
Beyond Leading Colour,” arXiv: 1505.01681 [hep-ph]. 

16. TLEP Design Study Working Group Collaboration, 
M. Bicer et al., “First Look at the Physics Case of 
TLEP,” JREP 1401 (2014) 164, arXiv: 1308.6176 
[hep-ex]. 

17. CMS Collaboration, S. Chatrchyan et al., “Observation 
of a new boson at a mass of 125 GeV with the CMS 
experiment at the LHC,” Phys.Lett. B716 (2012) 30-61, 
arXiv:1207.7235 [hep-ex]. 

18. ATLAS Collaboration, G. Aad et al, “Observation of a 
new particle in the search for the Standard Model Higgs 
boson with the ATLAS detector at the LHC,” Phys.Lett. 
B716 (2012) 1-29, arXiv: 1207.7214 [hep-ex]. 

19. ATLAS Collaboration, G. Aad et al., “Study of the spin 
and parity of the Higgs boson in diboson decays with 
the ATLAS detector,” arXiv: 1506.05669 [hep-ex]. 


20. CMS Collaboration, S. Chatrchyan et al., 

“Measurement of the properties of a Higgs boson in the 
four-lepton final state,” Phys.Rev. D89 no. 9, (2014) 
092007, arXiv: 1312.5353 [hep-ex]. 

21. A. Skjold and P. Osland, “Angular and energy 
correlations in Higgs decay,” Phys.Lett. B311 (1993) 
261—265, arXiv:hep-ph/9303294 [hep-ph]. 

22. B. Andersson, G. Gustafson, and B. Soderberg, “A 
Probability Measure on Parton and String States,” 
Nucl.Phys. B264 (1986) 29. 

23. T. Sjostrand, S. Ask, J. R. Christiansen, R. Corke, 

N. Desai, et al., “An Introduction to PYTHIA 8.2,” 
arXiv:1410.3012 [hep-ph]. 

24. C. Bierlich, G. Gustafson, L. Lonnblad, and A. Tarasov, 
“Effects of Overlapping Strings in pp Collisions,” JHEP 
1503 (2015) 148, arXiv:1412.6259 [hep-ph]. 

25. T. Biro, H. B. Nielsen, and J. Knoll, “Color Rope 
Model for Extreme Relativistic Heavy Ion Collisions,” 
Nucl.Phys. B245 (1984) 449-468. 

26. A. Bialas and W. Czyz, “Chromoelectric Flux Tubes 
and the Transverse Momentum Distribution in 
High-energy Nucleus-nucleus Collisions,” Phys.Rev. D31 
(1985) 198. 

27. DELPHI Collaboration, J. Abdallah et al, 
“Measurement of the Mass and Width of the W Boson 
in e+ e~ Collisions at </s = 161-GeV - 209-GeV,” 
Eur.Phys.J. C55 (2008) 1—38, arXiv:0803.2534 [hep-ex]. 

28. S. Catani, Y. L. Dokshitzer, M. Olsson, G. Turnock, 
and B. Webber, “New clustering algorithm for multi - 
jet cross-sections in e-b e- annihilation,” Phys.Lett. B269 
(1991) 432-438. 

29. ALEPH, DELPHI, L3, OPAL, LEP Electroweak 
Working Group Collaboration, J. Alcaraz et al, “A 
Combination of preliminary electroweak measurements 
and constraints on the standard model,” 

arXiv:hep-ex/0612034 [hep-ex]. 

30. A. Skjold and P. Osland, “Testing CP in the Bjorken 
process,” Nucl.Phys. B453 (1995) 3-16, 

arXiv:hep-ph/9502283 [hep-ph]. 

31. Y. Gao, A. V. Gritsan, Z. Guo, K. Melnikov, 

M. Schulze, et al, “Spin determination of 
single-produced resonances at hadron colliders,” 

Phys.Rev. D81 (2010) 075022, arXiv: 1001.3396 
[hep-ph]. 

32. S. Bolognesi, Y. Gao, A. V. Gritsan, K. Melnikov, 

M. Schulze, et al, “On the spin and parity of a 
single-produced resonance at the LHC,” Phys.Rev. D86 
(2012) 095031, arXiv:1208.4018 [hep-ph]. 

33. V. A. Khoze and T. Sjostrand, “Color correlations and 
multiplicities in top events,” Phys.Lett. B328 (1994) 
466—476, arXiv:hep-ph/9403394 [hep-ph]. 



